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Abstract. A recent revival of interest in synchronised oscillations due to neutrino-neutrino forward 
scattering in dense gases has led to two interesting applications with notable outcomes: (i) cosmo¬ 
logical bounds on neutrino-antineutrino asymmetries are improved owing to flavour equilibration 
prior to the onset of big bang nucleosynthesis, and (ii) a neutron-rich environment required for r- 
process nucleosynthesis is shown to be always maintained in a supernova hot bubble irrespective 
of flavour oscillations, contrary to results from previous studies. 1 present in this talk a pedagogical 
review of these works. 


1. THE NEUTRINO-NEUTRINO FORWARD SCATTERING SAGA 

It is well known that neutrino oscillations in a medium are affected by the presence of 
other particles. A familiar example is oscillations in the sun: Extra interaction channels 
available exclusively to the electron neutrino induce for it an excess “effective” mass 
through coherent forward scattering on the ambient electrons [1].^ Depending on the 
electron number density, this excess serves as a refractive index for the oscillating 
neutrino system, modifies the oscillation frequency and amplitude from their vacuum 
values, and is responsible for such phenomenon as the Mikheyev-Smimov-Wolfenstein 
(MSW) effect [3] and, with it, the large mixing angle (LMA) solution to the solar 
neutrino problem [4]. 

In certain astrophysical and cosmological settings, neutrinos form a gas so dense that 
forward scattering of neutrinos on the neutrinos themselves may constitute a significant 
source of refraction [5, 6]. Furthermore, such scattering differs markedly from scattering 
on other fermions, since in the former instance it is possible for the “beam” and “back¬ 
ground” neutrinos to exchange flavour via the neutral current. The net outcome is a set 
of “on-” and “off-diagonal” refractive indices that are dependent on the state of every 
other neutrino in the gas, thereby rendering the evolution of the whole ensemble highly 
nonlinear. 

In the very earliest works incorporating the refractive index from neutrino-neutrino 
forward scattering, the off-diagonal refractive indices were erroneously left out. This 
was rectified by Pantaleone in 1992 [7]. Shortly after, the implications of these “self in- 


^ This statement assumes that Ve oscillates into and/or v^. If a sterile neutrino is involved, then neutral 
current scattering of Ve on the surrounding electrons and nucleons will also contribute to the excess [2]. 
Note that we shall not be considering sterile neutrinos in this talk. 



teractions” for a multi-momentum gas were investigated extensively by Samuel [8], who 
diseovered in his numerieal studies that when a eertain eritieal density is reaehed, the gas 
begins to exhibit self-maintained eoherenee (a.k.a. synehronisation) when deeoherenee 
is a priori expeeted. This interesting effeet went unexplained for almost a deeade, and 
was reeently reexamined by Pastor, Raffelt, and Semikoz [9], who sueeeeded in giving 
it a elean and physieal interpretation on whieh I shall elaborate in this talk. 

The result of Ref. [9] has sinee been extended and applied to the only two environ¬ 
ments in whieh self interaetions are expeeted to play a dominant role: the early universe, 
and a eore eollapse supernova. In the first setting, synehronised oseillations are investi¬ 
gated in the eontext of relie neutrino-antineutrino asymmetries in the epoeh preeeding 
big bang nueleosynthesis (BBN). Independent studies by several groups have demon¬ 
strated that if the neutrino oseillation parameters are indeed those inferred from the 
solar and atmospherie neutrino data, then the asymmetries in the and the seetors 
must be equilibrated with that in the Vg seetor prior to BBN, and be subjeet to the same 
stringent constraints on the Vg asymmetry derived from uncertainties in the primordial 
^He abundance [10, 11, 12] (see also Ref. [13]). In the second case, self interactions 
modify in an unexpected way the flavour evolution of both neutrinos and antineutrinos 
in the hot bubble region above a nascent neutron star a few seconds after core bounce 
[14]. This impinges directly on the success or otherwise of r-process nucleosynthesis, 
which is believed to take place in this setting. 

In this talk, I shall present a pedagogical review of these works, starting with the inter¬ 
pretation of synchronised oscillations in Sec. 1, and then moving on to the applications 
in Secs. 2 and 3. The material used here is drawn largely from Refs. [9], [10], [11], [12], 
and [14]. If you the reader happen to the author of one of these works, please pardon me 
for quoting you almost verbatim in some places. 


1.1. Equations of motion 

1.1.1. Density matrix formalism 

We begin with a simple two-flavour system. In order to understand the synchronisation 
phenomenon, it is convenient to work with bilinears of the neutrino wave function (and 
hence density matrices), since, as we shall see, the neutrino self interaction potential 
appears also in this form. 

For a group of neutrinos with momentum p, the density matrix at time t in the flavour 
basis is defined as 
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where the subscripts a and j5 denote two different neutrino flavours, Nq{p) is some 
time-independent normalisation factor, and 

\¥i{P,t)) = ai{p,t)\Va) +bi{p,t)\Vp) ( 2 ) 

is the wave function of the ith neutrino in the group. The functions |a-(p,t)p and 
\b-{p^ t) p are probabilities in the one-particle formulation. Thus the sums 

Nv^=Y^\ai{P:t)\^ = PaaNQ{p), Ny^=Y,\b^{p,t)\^ = p^^NQ{p) (3) 

i i 


in the diagonal of p{p,t) represent respectively the occupation numbers of and 
at a given momentum, while the off-diagonal p^^ and are intangible quantities 
containing the system’s phase information. 

We evolve p(p,t) by taking its commutator with the Hamiltonian H{pp) in flavour 
space, 

dtp{pp) = -i[H{p,t),p{p,t)]. (4) 


As an example, the Hamiltonian governing Vg .j; oscillations in the sun takes the 
form 
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in which the first term describes vacuum oscillations, and the second term with matter 
potential A (p,t) = \/lGpn^_ (t), where {t) denotes the solar electron number density, 
is responsible for refractive effects.^ 

For our purposes, it is useful to recast the equation of motion (4) into Bloch form by 
first parameterising p{p,t) and H{p,t) in terms of the functions (/q, P) and (Vq, V): 
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Here, o = OxX + ovy -f cr^z are the Pauli matrices. The quantities P = PxX + Pyy + PzZ 
and V = \4x -|- Vyy -|- VzZ may be regarded, respectively, as a three-dimensional “spin 
polarisation” and a “magnetic field” vector, whose time and momentum dependences 
are implicit. The equation of motion (4) transforms into 


<9,P = VxP, o>,Po = 0, (7) 

such that the dynamics of the oscillating system is analogous to a spin vector predisposed 
to processing around a magnetic field V at a rate |V|, with 


Nvg,iPP)-Ny^ip,t) 


( 8 ) 


^ Rigorous derivations of the equation of motion and various interaction potentials can be found in Refs. 
[15] and [16]. See also Ref. [17]. 



parameterising the excess of over v^, and conservation of probability demands that 
the function 




(9) 


remain constant with time. 

Observe in Eq. (7) that Vq is redundant for flavour oscillations, since it contributes only 
an unobservable overall phase to the neutrino wave function. Some typical examples of 
the field vector V are (i) 
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for oscillations in vacuum, where, pictorially, the unit vector B forms an angle 29 with 
the negative z-axis and 5m^/2p is the precession frequency, and (ii) 
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for oscillations in a medium, with the matter potential A(p, t) entering in the z-direction. 
As mentioned before, the presence of A(p, t) generally results in the modification of the 
oscillation frequency and amplitude. In the spin precession picture, this is equivalent to 
defining a field vector B^^^ 
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which makes an angle 29^^ with the negative z-axis, and P precesses around is at a rate 
5m^^^j2p, where 


sin20^ff 
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{2pA{p,t) / 5m^ — cos29)^ + sin^ 29 
5n^\J {2pA{p, t) / dm^ — cos 20)^ + sin^ 20. 


(13) 


If \A{p,t)\ is much larger than |5m^/2p|, the former will force the total field vector 
to align with the z-axis. The effective mixing angle ^eff approaches 0 or ;r/2, thereby 
suppressing flavour oscillations. On the other hand, if the condition 
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happens to hold, then the vector V points in the x-direction, forming a right angle 
with the z-axis which corresponds to maximal mixing, 0^^^ = Tt/A. This is the so-called 
resonance condition. 

Finally, observe that the field vector B^^^ differs for each momentum mode. Thus the 
spin vectors of an ensemble of neutrinos possessing a spectrum of momenta do not in 
general precess in the same direction, or at the same rate. 



1.1.2. Add some neutrino-neutrino forward scattering 


Suppose now we have a spatially homogeneous and isotropie neutrino gas of volume 
V. Contribution from neutrino-neutrino forward seattering to the total field veetor V 
appears in the form [16] 


V2Gf- 
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J, 


(15) 


where 

J=EPA(Pi) (1® 

k 


is the total polarisation veetor, and the subseript k labels the kth momentum mode. Thus 
the simplest equation of motion one ean write down to deseribe flavour oseillations in 
the kth momentum mode of a self-interaeting neutrino gas will eonsist of a term from 
vaeuum oseillations, and another from self interaetions: 


d,Pt = 
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Note that Eq. (17) is valid only for the ease of aetive-aetive flavour oseillations. For 
an aetive-sterile system, all interaetion potentials must eouple exelusively to the aetive 
neutrino, and their eontributions to the total field veetor V in the Bloeh formulation 
always point in the ±z-direetion. 


1.2. Physical interpretation of synchronised oscillations 

1.2.1. Vacuum-like oscillations 


In the pioneering work of Samuel [8], Eq. (17) was solved numerieally for a neutrino 
gas with a distribution of momenta. To his surprise, the oseillating multi-momentum gas 
exhibited synehronised, monoehromatie behaviours when the self interaetion potential 
was allowed to dominate over the vaeuum oseillation rate: 


V^Gp,^, dnd 

^ ^Pk 


(18) 


A series of papers by Samuel and eo-workers on the same subjeet emerged in the 
years following this initial diseovery [18, 19, 20], extending the results of Ref. [8] to 
inelude antineutrinos and a non-neutrino baekground in the eontext of the early universe. 
Exaet solutions to the equations of motion were found [19], and the stability of the 
synehronised state analysed [20]. Nonetheless, a elean, physieal explanation for the 
synehronisation phenomenon had eontinued to elude the players. This interesting effeet 
was elueidated reeently by Pastor, Raffelt, and Semikoz in Ref. [9], whose argument 
goes as follows. 



When the neutrino gas is suffieiently dense and the eondition (18) satisfied, the 
equation of motion (17) for a single momentum mode becomes dominated on the right 
hand side by the self interaction potential: 


P) P 


This implies an approximate solution 
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( 19 ) 


( 20 ) 


where the superscript “i” denotes initial. In plain English, all polarisation vectors P^ 
across the momentum distribution are individually first and foremost drawn towards the 
giant total polarisation vector J, and process around it at a rapid rate of y2G^|j|/y.if 
the evolution of J proceeds at a rate much slower than this precession frequency, then all 
P^’s will remain pinned to J, which, as we shall see, forms the basis of synchronisation. 

The dynamics of J is governed by an equation of motion constructed from summing 
the single mode equation (17) over all momentum modes. 


k ^Pk 


( 21 ) 


Substituting the approximate solution (20) and averaging the oscillations to zero over 
the momentum spread, we immediately see that the evolution of J is non-dissipative. 
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\^k^k^Q^Pk)\ 

Since the slowly-moving part of every polarisation vector P^ follows to first approxi¬ 
mation the dynamics of J, an apparently monochromatic behaviour emerges within the 

multi-momentum neutrino gas, and is the synchronised oscillation frequency 

with an explicit dependence on the initial conditions. 


7 . 2 . 2 . Add a non-neutrino background 


Granted that neutrino-neutrino forward scattering continues to constitute the domi¬ 
nant interaction potential, it is easy to generalise the previous picture to include a non¬ 
neutrino background A(p,t). The equation of motion is essentially the same as Eq. (17), 
save for the replacement 
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and the analogue of the “dense gas” condition (18) is 
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When this condition is satisfied, all polarisation vectors become pinned to J as per 
the approximate solution (20), and the evolution of J is equally well described by Eq. 
(22) upon replacing 


5m^ \ / dm^ 

-iFr^KW 
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where 






(27) 


represents an average background. 

Perhaps the most interesting aspect of this result is that a strong self interaction 
potential not only forces all modes in a multi-momentum neutrino gas to oscillate at the 
same frequency, but also with the same amplitude, regardless of the individual modes’ 
response to the non-neutrino background. Previously, we saw that electrons in the solar 
interior alter the mixing angle and the oscillation frequency for Vg oscillations in 

a momentum-dependent fashion. In the present case, however, the effective field vector 
for J in Eq. (26) necessarily implies a common mixing angle 


sin 20 


synch 
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\/((2l(p,O)/(^)-cos20)2 + sin220 


(28) 


and a common oscillation frequency 


"synch = (^) y^((^(P,O)/^^^-cos20^ +sin220, (29) 

for all momentum modes. A corollary of Eqs. (28) and (29) is that if the “average” 
resonance condition 

{Mp)) = ( y cos20 (30) 

is met, then it is necessarily met by all momenta at the same time or at the same 
place, and any MSW resonant conversion that is incurred must proceed with the same 
adiabaticity across the spectrum [11, 12]. 


1.2.3. Add some antineutrinos 

In a realistic system such as the early universe, it is inevitable that neutrinos and an¬ 
tineutrinos will coexist, in which case it is necessary to consider also neutrinos scattering 



on antineutrinos and vice versa. To this end, we first define the density matrix for an¬ 
tineutrinos of momentum p in flavour space, and the corresponding polarisation vector 

P: 
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To first order in the equations of motion are [16] 


x/^Cj — 

— — — 
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where J is the antineutrino analogue of J, and we have segregated the background 
interaction potentials into a CP symmetric (A^^^^) and a CP asymmetric (A^p ) part, 
with implicit momentum and time dependences. A CP asymmetric potential, such as 
that arising from a local four-fermion coupling, carries opposite signs in the interaction 
Hamiltonian for neutrinos and antineutrinos, while is identical for both v and v, 
and comes usually from higher order propagator effects [5, 16]. 

Synchronisation occurs when the condition 
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holds, and the vector I = J — J now assumes the role played by the total polarisation 
vector J in the previous neutrino-only case. All P^ and P^ become pinned to I, 


P^ ~ (Pj,-I^)I-f precessions around I, 

P^ ~ (P^-1^)1-I- precessions around I, 

and the corresponding evolution equation for the compound system is 


where 
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Observe that the present formulation is invalid if 

k q 

although the ease of identieal v and v speetra ean lead to some eurious effeets [9]. 
It is interesting to note that both neutrinos and antineutrinos are synehronised and 
evolve in an identieal manner, regardless of the presenee of a baekground medium that 
distinguishes between the CP partners [11, 14]. This has an interesting implieation for 
r-proeess nueleosynthesis in a supernova hot bubble, to be diseussed later. 

This is as mueh baekground material as we shall need. Let us now turn to the 
applieations. 


2. APPLICATION 1: EQUILIBRATION OF RELIC 
NEUTRINO-ANTINEUTRINO ASYMMETRIES 

One of the many open questions in eosmology is the possibility of admitting a large relie 
neutrino-antineutrino asymmetry. 




— Hr 


(38) 


where ny, = f Ny,p^dp denotes the number density of the partiele speeies t/r. At temper¬ 
atures T > 2 MeV, thermal and ehemical equilibria are generally believed to hold for 
an ensemble of neutrinos and antineutrinos of some aetive flavour a. The distribution 
funetions and A,-, assume the Fermi-Dirae form, 
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where 
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and ~ asymmetry in v« and rnay be alternatively expressed 

in terms of the speeies’ ehemieal potential (^y^. 
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with Uy = 2^(3)r^/and ^ is the Riemann zeta funetion.^ 


^ The terms neutrino-antineutrino asymmetry, neutrino chemical potential, and neutrino degeneracy 
parameter will be used interchangeably throughout this talk. 



As of now, there are no direct observations of the cosmic neutrino background. Thus 
the existence or otherwise of any sizeable ^ can only be established indirectly from the 
study of the cosmic microwave background radiation (CMBR) anisotropy spectrum, and 
from requiring consistency with the observationally highly successful theory of BBN. 
In the former case, we note that the contribution by one neutrino species (v« and Va) to 
the total energy density in radiation in the universe depends immediately on the species’ 
chemical potential: 


V„v„ 






(42) 


An increase in the energy density in radiation tends to delay the epoch of matter- 
radiation equality, which is known to magnify the amplitude of the first acoustic peak in 
the CMBR angular power spectrum, and shift all other peaks to higher multipoles [21]. 
Currently, CMBR limits ^ to below 3, applicable to all neutrino flavours [22]. 

Secondly, an asymmetry in the electron neutrino sector impinges directly on the j8- 
processes 

n-\-Ve ^ p + , p + Ve^n + e^ (43) 

that are responsible for setting the neutron-to-proton ratio prior to BBN (T ~ 1 MeV)— 
a positive tends to lower n^/np while a negative ^y^ raises it. Since the primordial 
^He mass fraction Yp depends exclusively on the value of nnjnp at the weak freeze-out, 
the measurement of Yp can be used to limit the allowed values of ^y^, which cover the 
range —0.01 < < 0 .07 at present [10, 23]. Furthermore, a nonzero ^ in any neutrino 

flavour tends to accelerate the expansion of the universe through its contribution to the 
energy density a la Eq. (42). This has the effect of raising the freeze-out temperature 
for the processes of Eq. (43) and therefore also the instantaneous neutron-to-proton 
ratio, affecting again the value of Yp. An interesting nucleosynthesis scenario that allows 
for large neutrino asymmetries while simultaneously preserving the standard primordial 
^He yield is known as degenerate BBN (DBBN), whereby an increased due to a 
large ^y^ and/or ^y_^ is compensated for by a sizeable positive ^y^ [23]. Several recent 
combined analyses of (D)BBN and CMBR have generated the constraints 

-0.01 < < 0.22, |^v„vj<2.6, (44) 

assuming no neutrino oscillations [24, 25]. 


2.1. Neutrino oscillations? 

Naturally, one is curious to find out how the numbers in Eq. (44) hold up in the pres¬ 
ence of oscillations. To date, there are three pieces of evidence for neutrino oscillations, 
(i) Electron neutrinos originating from nuclear reactions in the sun have been observed 
to transform largely into and/or [26]. Currently, these oscillation parameters all 
provide acceptable explanations for the Vg deficit [4]: 


EMA 


~ 5 X 10 ^ eV^ 


sin20 ~ 0.8, 



EOW 

5n? ~ 10 ^ eV^ 

sin 20 ~ large, 


Vacuum 

5n? ~ 10^ eV^ 

sin 20 ~ large, 


SMA 

8rr? c:=l1 X 10 ^^ 

sin20 ~ 0.05. 

(45) 


The LMA solution, however, is unanimously the most favoured, (ii) The disappearance 
of atmospheric and is generally attributed to oscillations into and v^, with oscil¬ 
lation parameters 5m^ ~ 3 x 10~^ eV^ and sin20 ~ 1 [27]. (iii) The Liquid Scintillator 
Neutrino Detector (LSND) experiment has reported evidence for Vg oscillations 
with 5m^ > 1 eV^ and sin 20 ~ 0.03 —0.1 [28]. Since it is impossible to accommodate 
all three signals with only three neutrinos, the most common practice is to either intro¬ 
duce a fourth sterile neutrino (which in fact does not provide a very good fit to the data 
[29]), or ignore the LSND result. We shall be considering the latter option for simplicity. 

In a rudimentary study by Lunardini and Smirnov [13], it was suggested that three 
flavour oscillations in a near-bimaximal mixing framework (i.e., if LMA holds) nec¬ 
essarily lead to flavour equilibrium prior to BBN (see also an earlier work by Savage, 
Malaney, and Fuller [30]). This proposal was later examined in detail in the numerical 
studies of Dolgov et al. [10], and then confirmed by analytical means in a paper by my¬ 
self [11], and in another independent work by Abazajian, Beacom, and Bell [12]. The 
common conclusion of these works is that constraints from BBN on the Vg asymmetry 
must now apply also to asymmetries in the and the sectors, which constitutes a 
significant improvement over the DBBN/CMBR bounds quoted in Eq. (44). 

A rigorous three flavour analytical study is at present beyond our means. Fortunately, 
the hierarchical nature of the mass splittings allows us to understand the full equilibration 
mechanism in terms of two separate equilibrating transitions—the first between Ly 
and Ly^, governed by the atmospheric neutrino oscillation parameters, and the second 
between Ly^ and Ly^, where denotes some linear combination of and and the 
Ve Vjc transformation proceeds with the parameters of Eq. (45).^ I shall outline the 
mechanics of these transformations. 


2.2. Equilibration of Ly^ with Ly^ 

Suppose there is some pre-existing asymmetry in either or both of the Vg and the Vy 
sectors. The one-body reduced density matrices for neutrinos and antineutrinos are as 
defined previously. These are again parameterised in terms of two polarisation vectors 
P and P, and for the present application, we choose the normalisation factor Nq{p) to 
be A^eq(O) as defined in Eq. (40). The evolution of the neutrino-antineutrino ensemble 
at temperatures T ~ ^(1) —i^(lO) MeV is governed by Eq. (32), where the volume V 
accompanying the self interaction term is already absorbed into the present definition of 


^ In principle, if the mixing angle commonly known as Pjj is nonzero, a third equilibrating transition is 
also present. I shall not talk about this possibility here. The interested reader is referred to the original 
works for details. 



Nq{p), and the background potentials are 


^CP+ ~ 


3ni^ 


■^eei 


^QP_ — y/^GpnyP, 


(46) 


in which Eee is the electron-positron energy density, m^y is the mass of the W boson, 
and 77 is the electron-positron asymmetry. By universal charge neutrality, 77 must be of 
the order of the baryon asymmetry so that is negligible in comparison with 
The quantity we are interested in is the difference between the two asymmetries 
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Hence we would like to track the evolution of the vector I = (J — J)/ny (note the new 
definition), and flavour equilibrium necessarily implies 4 = 0 - 

The next question is, are all the modes synchronised so that we can apply results 
from the previous sections? The answer comes from figuring out the conditions under 
which Eq. (33) holds true, and for the temperature range of interest, this amounts to 
requiring only that the disparity between the initial asymmetries be larger than, say, 
10 ^ in magnitude [11]. Then, assuming that all neutrinos and antineutrinos are initially 
in thermal and chemical equilibrium, i.e.. 


pi ti ^ ^eq(‘^vj 77eq(<^y^) _ ■ ^ A/eq( 

fVeq(O) ’ ■ ■ iVeq(O) ’ ^ ^ 

the general compound evolution equation (35) can be rewritten explicitly for this appli¬ 
cation as 


a,I 
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SV2G^p ] 
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xl, 


(49) 


where 



(50) 


represents some average momentum. 

Equation (49) has a straightforward interpretation. Consider the terms inside the 
parentheses. These are identically the vacuum and electron-positron background terms 
that one would find in a single momentum evolution equation, and control the ensemble’s 
synchronised mixing angle. 




sin 20 


{IpAppp/ 5m2 — cos 20 ) 2 -|- sin^ 20 


( 51 ) 


p=p 



evaluated for = —SV2GppEee/3m^ at p = p. The faetor multiplying the braek- 
eted terms, whieh we shall label as k, 

does not play a role in determining the effeetive mixing angle. It does, however, influenee 
direetly the oseillation frequeney. 


^»^synch 

2p 


Kdtrr\ {IpApp,/8m?- — eos20)2 + sin^20 
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(53) 


For instanee, vanishes for = —<^v^, and oseillations are switehed off eom- 

pletely. 

An approximate solution to Eq. (49) ean be obtained in the adiabatie limit. Speeifi- 
eally, the expression for the variable 4 as a funetion of time is 


assuming the validity of the adiabatie eondition 


VzdtVx-VxdtVz 

fC(y2 + y^2)3/2 


(55) 


where 14 = {8n? jlp') sin20, and V), —?^cp+ ~ t^P) cos 20, both evaluated for p = 
p. Equation (54) prediets for maximal vaeuum mixing an MSW-like effeet, transforming 
4 from 4 to 0 (plus some small amplitude oseillations) when vaeuum oseillations 
overcome refraetive matter effeets. The temperature at whieh this equilibrating transition 
takes plaee ean be established roughly by solving 


(56) 

3m^ Ip 

Eor initial ehemieal potentials satisfying the eonstraints (44), the equilibration tempera¬ 
ture turns out to be ~ 2.6 MeV for the EMA solution. In the eases of the EOW and the 
Vaeuum mass splittings, the temperatures are ~ 0.9 MeV and ~ 0.3 MeV respeetively. 
Evidently, only the EMA equilibrating transition ean take plaee well ahead of BBN.^ 
Eor the oseillation parameters of the SMA solution, 4 remains elose to its initial value 
even after the “transition” at T ~ 1.9 MeV, sinee both the vaeuum and the baekground 


^ The real LMA solution encompasses a range of mixing parameters that are merely large, but not 
maximal [4]. Thus only a partial equilibrium between the asymmetries can be achieved from synchronised 
oscillations alone. However, collisions with the background medium (i.e., momentum-changing non¬ 
forward scattering) are inevitable, and these are expected to drive the equilibration to a more complete 
state. 



terms are predominantly in the negative z-direction for dnp- > 0, and the usual MSW 
resonance condition cannot be satisfied. 

A second deciding factor on the efficacy of flavour equilibration is the adiabaticity of 
the transition from matter-suppressed to vacuum oscillations. Unlike that encountered 
in, for instance, solar neutrino analyses, the adiabaticity parameter 7 of Eq. (55) is 
strongly dependent on the initial conditions. Take for concreteness the case of = 
0. The adiabatic condition 7 < 1 always holds for LMA if > 0.01, but can be 
badly violated at the transition point for the LOW 5m^ unless > 0.1 [11]. In 
the case of maximal mixing, violation of adiabaticity at the transition point generally 
results in large amplitude “post-transition” oscillations about the equilibrium point at 
an angular frequency roughly equal to 5»^synch/^^- Naturally, this is quite a separate 
phenomenon from true equilibration. On the other hand, an adiabaticity parameter that 
evaluates to infinity at all times (because, for example, k ^ 0) signifies that there is no 
transition at all. From the perspective of equilibrating two vastly different asymmetries, 
the requirement of \^yj > 0.01 (assuming = 0) in the LMA case for a smooth 
transition is, by definition, not a major concern. 


2.3. Equilibration of Ly^ and Ly^, and new constraints on 
neutrino-antineutrino asymmetries 

Super-Kamiokande’s atmospheric neutrino data indicate maximal mixing between 
and V^:, with a mass splitting that is some two orders of magnitude larger than the 
maximum acceptable solar 5m^. This means, in general terms, that equilibration of the 
asymmetries in the and the sectors will occur in a manner similar to that outlined 
before, but at a higher temperature, and subject to a background potential proportional to 
the muon-antimuon energy density . A comparison of the magnitudes of the vacuum 
and the background potential terms returns an equilibration temperature of ~ 12 MeV 
for 5m^ = 3 X 10^^ eV^ and p = p k/T, and complete equilibrium between Ly^ and 
Ly^ can always be achieved because of maximal mixing and collision-induced flavour 
relaxation not discussed here [ 10 , 11 , 12 ]. 

Thus, summing up the findings, we conclude that all active neutrino-antineutrino 
asymmetries or chemical potentials must come to an equal value before the onset 
of BBN, if the flavour oscillation parameters are those indicated by the atmospheric 
neutrino data and by the solar LMA solution. The stringent constraint imposed by the 
primordial ^He abundance on ^y^, 

?v, - Iv, - iv, < 0.07, (57) 


now applies to all three neutrino flavours. 



3. APPLICATION 2: /^-PROCESS NUCLEOSYNTHESIS IN A 

SUPERNOVA HOT BUBBLE 


Heavy nuclei with mass number A > 70 are predominantly produced by slow and 
rapid neutron capture—the s- and the r-processes. The most plausible site for r-process 
nucleosynthesis suggested so far is the hot bubble between a protoneutron star and the 
escaping shock wave in a core collapse supernova a few seconds after core bounce. A 
key condition for this process to be successful, however, is that the environment must 
be rich in neutrons (relative to protons), and neutrinos released from the cooling of the 
nascent neutron star play an important role in establishing this condition. 

Just as in the early universe, the neutron-to-proton ratio is fixed by the -processes in 
Eq. (43). Assuming charge neutrality, it is common to express this ratio in terms of the 
electron fraction Yg'. 


Hi - 1. 

Yg 


(58) 


Clearly, Yg < 0.5 will satisfy the minimum requirement for neutron predominance, but 
a successful r-process may call for Yg < 0.45. Because neutrinos are far more abundant 
than electrons and positrons, Yg is governed largely by the neutrino spectra and fluxes. 
Near the weak freeze-out radius (~ 30 —35 km), this is well approximated by [31] 





(59) 


where Ly^ and are the electron neutrino and antineutrino luminosities (typically of 

order 10^^ erg s"^), £ = {ElJ/{Ey^), and £ = (E? )/{Ey ). Neutrinos streaming out of 
the neutrino sphere at a radius of ~ 10 km generally possess spectra of Fermi-Dirac 
form, with mean energies 


{EyJ ~ 11 MeV, {E^J ~ 16MeV, 


{E^ ^ V V ) - 25 MeV. 


(60) 


These numbers yield for Yg a value of ~ 0.41, allowing for a successful r-process, in the 
absence of neutrino oscillations. 

Contrastingly, if flavour oscillations between Vg and occur outside the neutrino 
sphere but still within the weak freeze-out radius, then the two dissimilar neutrino 
spectra will be exchanged to various degrees depending on the oscillation parameters, 
and Yg becomes modified as a generic consequence [31]. Since a large, positive electron- 
positron asymmetry is always present in the background medium, and the interaction 
Hamiltonian dominated by the CP asymmetric potential \/2Gp {n^_ — ), conventional 

wisdom tells us that the MSW resonance condition cannot be satisfied simultaneously by 
both neutrinos and antineutrinos. If the former meets the condition, the more energetic 
V^’s will be transformed across the resonance into Vg’s, while the Vg energy spectrum 
remains virtually unaffected. This then leads to an over-production of electrons, pushes 
Yg to above 0.5, and consequently ruins the r-process. The mass splitting required for 
such oscillations to be effective is 5m^ > leV^, which overlaps with the region of 
parameter space needed to explain the LSND result [28]. 



This situation changes drastically when neutrino-neutrino forward scattering is also 
taken into account. Specifically, for the nominal neutrino luminosity Ly = 10^^ erg, 
the electron fraction always stays below 0.5, such that the minimum requirement for 
a neutron-rich environment can always be satisfied, irrespective of neutrino oscillations. 
This interesting result can be easily understood in terms of synchronisation [14]. 

As noted before, when neutrino-neutrino forward scattering constitutes the dominant 
interaction potential, it synchronises both neutrinos and antineutrinos, regardless of the 
nature of the non-neutrino background medium. Consequently, if the neutrinos were to 
encounter an MSW resonance, the antineutrinos would also be dragged along for the 
ride, and any spectral swap that accompanies the transitions would occur with identical 
efficiency. The net effect is that the post-resonance Vg and Vg will now have identical 
spectra, having taken over those which once belonged to and respectively. 
Assuming that Ly^ ~ Ly^, Eq. (59) always gives Yg ~ 0.5 for the electron fraction. 

4. CONCLUSIONS 

I have reported in this talk the findings of a series of recent works devoted to explaining 
and applying the phenomenon of synchronised flavour oscillations due to neutrino self 
interaction, in the early universe and in a core collapse supernova. 

For the “explanation” part, I presented the original interpretation of Pastor, Raffelt, 
and Semikoz [9], as well as extensions to it introduced by myself [11], and separately, 
by Abazajian, Beacom, and Bell [12]. 

Application-wise, the analyses of Dolgov et al, and Refs. [11] and [12] demonstrate 
that neutrino self interactions in the epoch immediately prior to BBN are sufficiently 
intense to synchronise the multi-momentum flavour oscillations, which, together with 
the large mixing angles indicated by present neutrino data, allows for the efficient 
equilibration of all pre-existing neutrino-antineutrino asymmetries. I have reviewed 
these works in considerable detail, and relayed the unanimous conclusion that stringent 
BBN bounds on the electron neutrino-antineutrino asymmetry must now apply also to 
asymmetries in the and the V^: sectors. 

Lastly, I spoke briefly about a possible upset of the neutron-rich environment and thus 
the success of r-process nucleosynthesis in a supernova hot bubble by “conventional” 
neutrino oscillations. I proceeded to report the new results of Pastor and Raffelt [14], 
which show that, if neutrino self interactions are properly taken into account, synchro¬ 
nised oscillations will always prevent the neutron-to-proton ratio from dropping below 
the minimum value required for neutron predominance. 

And on this note, I rest my case. 
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